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Abstract 

In this paper we study the possibly existing anisotropy in the accelerating expansion Universe 
by use of the full sample of Union2 data. Using the hemisphere comparison method to search 
for a preferred direction, we take the deceleration parameter qo as the diagnostic to quantify the 
anisotropy level in the wCDM model. We find that the maximum accelerating expansion direction 
is (Z,6) = (314°~2oo )i with the maximum anisotropy level of Aqo,max/qo = 0.79^0 27, and 

that the anisotropy is more prominent when only low redshift data {z < 0.2) are used. We also 
discuss this issue in the CPL parameterized model, showing a similar result. 



PACS numbers: 98.80.Es, 98.80.Jk 

1. INTRODUCTION 

Since the discovery of cosmic acceleration in the late 
of last century |1], the type la supernovae (SNIa) have 
become an important tool in determining the cosmolog- 
ical parameters. Given certain cosmological models, one 
usually implements the joint analysis with SNIa in combi- 
nation with other observations, such as large scale struc- 
ture 12], the cosmic microwave background (CMB) radi- 
ation [5], and so on, to give constraints of cosmological 
parameters. 

Two pillars of modern cosmology are general rela- 
tivity and the cosmological principle. It is assumed 
that Einstein's general relativity still holds on the cos- 
mic scale, which means that the evolution of the uni- 
verse is governed by general relativity. The cosmological 
principle [5j says that our universe is homogeneous and 
isotropic on the cosmic scale. Indeed, the assumption 
of homogeneity and isotropy is consistent with currently 
accurate data coming from the cosmic microwave back- 
ground (CMB) radiation, especially from the Wilkinson 
Microwave Anisotropy Probe (WMAP) [S], the statistics 
of galaxies [6], and the halo power spectrum [7j, etc. And 
current astronomical observations are in agreement with 
AC DM model [Sj. 

However, the standard model is also challenged by 
some observations [9] (for more details see [10^ and refer- 
ences therein) . Thus it is necessary to revisit the assump- 
tion of the homogeneity and isotropy. As more and more 
supernovae data are released [HI [12] , this study becomes 
possible. 

From the theoretical point of view, the anisotropy 
may arise in some cosmological models. For example. 



a vector field may lead to anisotropy of the universe and 
gives rise to an anisotropic equation of state of dark en- 
ergy [13|. Peculiar velocities are also associated with 
dipole-like anisotropics, triggered by the fact that they 
introduce a preferred spatial direction, and as a result 
of the drift motion, observers may find that the accel- 
eration is maximized in one direction and minimized in 
the opposite [Hj. By use of 288 SNIa [H], Davis et al. 
studied the effects of peculiar velocities on cosmological 
parameters, including our own peculiar motion, super- 
nova's motion and coherent bulk motion, and found that 
neglecting coherent velocities in the current sample would 
cause a systematic shift in the equation of state of dark 
energy, with deviation Aw = 0.02 [16] . In this paper, 
we focus on the issue related to the possible existence of 
anisotropy, and search for a preferred direction by using 
the SNIa Union2 (consisting of 557 SNIa) data^. Some 
previous works payed attention to this issue using the 
SNIa data and found no statistically significant evidence 
for anisotropics 17 . Using the Union2, some authors de- 
rived the angular covariance function of the standard can- 
dle magnitude fluctuations, searching for angular scales 
where the covariance function deviates from in a statis- 
tically significant manner, and no such angular scale was 
found [18] . Yet on the other hand, using the SNIa data in 
the frame work of an anisotropy Bianchi type I cosmolog- 
ical model and in the presence of a dark energy fluid with 
anisotropy equation of state, Ref. jTHI found that a large 
level of anisotropy was allowed both in the geometry of 
the universe and in the equation of state of dark energy. 
Ref. [10] constructed a "residual" statistic, which is sen- 
sitive to systematic shifts in the SNIa brightness, and 
used this to search in different slices of redshift for a pre- 
ferred direction on the sky, and found that at low redshift 
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{z < 0.5) an isotropic model was barely consistent with 
the SNIa data at 2-3a. In addition, Ref. [21] took use 
of the hemisphere comparison method to fit the ACDM 
model to the supernovae data on several pairs of opposite 
hemispheres, and a statistically significant preferred axis 
was found. More recently, Antoniou and Perivolaropou- 
los f22J have applied the hemisphere comparison method 
to the standard ACDM model and found that the hemi- 
sphere of maximum accelerating expansion is in the di- 
rection {l,b) = (309°+|3o , 18°+i°o° ) with Union2 data. 
This result is consistent with other observations, such 
as CMB dipole [13], CMB quadrupole [H], CMB oc- 
topole [3S], large scale velocity flows [2S] and large scale 
alignment in the QSO optical polarization data [27]. In 
these observations, all the preferred directions appear 
to be towards the North Galactic Hemisphere. They 
obtained the average direction of the preferred axes as 
{l,b) = (278° ±26°, 45° ±27°). 

In this paper we study the dependence of this result 
on the dark energy model. We consider two models. One 
is the wCDM model; the other is a dynamical dark en- 
ergy model with the CPL parametrization [55]. We take 
the present value of deceleration parameter qo as the di- 
agnostic to quantify the anisotropy level of two opposite 
hemispheres. 

The paper is organized as follows. In the next section 
we give a general introduction to the hemisphere com- 
parison method using qo to quantify the anisotropy level 
and we apply it to wCDM model fitted by the Union2 
dataset. In Sec. 3, we give the numerical results on the 
preferred directions from the SNIa data with different 
slices of redshift. We also give the result for the CPL 
model. Sec. 4 gives our conclusions. 



HEMISPHERE COMPARISON METHOD 
USING THE UNION2 DATASET 



Since its first release, the SNIa data have been becom- 
ing an important tool for us to understand the evolution 
of the universe. And as time goes on, more accurately de- 
termined data will be released. For example, the future 
Joint Dark Energy Mission ( JDEM) ^ is aimed to explore 
the properties of dark energy and to measure how cos- 
mic expansion has changed over time. In this paper we 
take use of the Union2 dataset [IT, which contains 557 
type la SNIa data and uses SALT2 for SNIa light-curve 
fitting, covering the redshift range z = [0.015,1.4] and 
including samples from other surveys, such as CfA3 ^29j. 
SDSS-II Supernova Search |30j and high-z Hubble Space 
Telescope. 

We fit the SNIa data by minimizing the value of 
the distance modulus. The xin foi' SNIa is obtained 
by comparing theoretical distance modulus ^th{z) = 



h\ogio[dL{z)] ± /Lto, where = 42.384 — blog^oh is a 
nuisance parameter, with observed /^o& of supernovae: 



i=l 



For a flat FRW cosmological model, one has 



In the case of wCDM model, we have 

H\z) = H^[n„M'^ + zf + {i-n„,o){i + zf+^n- (2.2) 

And the present value of the deceleration parameter can 
be expressed as 



1 3 



(2.3) 



In the case of the CPL parametrization, the equation of 
state of dark energy is w ^ wo + wi . Accordingly one 
obtains 

H{z) - H^[n^o{l + zf + {l-n„,o){l + zf+^'^»+'^' 
exp(-3wiz/(l + z))], (2.4) 
1 3 



-wo(i - 



7n0) 



(2.5) 



Since /ig is a nuisance parameter, we can eliminate the 
effect of ^0 in the following way. We expand xin with 
respect to fio [ST] : 



A + 2Bfj,o + Ci4, 



(2.6) 



where 



A 
B 
C 
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Eq. (2.6l has a minimum as 



Xsn 



X 



sn.min 



A - pyc, 
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which is independent of /io. In fact, it is equivalent to 
performing a uniform marginalization over /ig, the dif- 
ference between and the marginalized xin j^^t a 
constant [3T]. We will adopt xln the goodness of fit- 
ting between theoretical model and SNIa data. 

The directions to the SNIa we use here are given in 
Ref. [T5] , and are described in the equatorial coordinates 
(right ascension and declination). In order to use the 
hemisphere comparison method, we need to convert these 
coordinates to usual spherical coordinates (0,0), and fi- 
nally to the galactic coordinates {l,b) ,32]. 
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We apply the hemisphere comparison method to find 
the possibly existing preferred axis. This method was 
first proposed in Rcf. [2T], and further developed in 
Ref. '52 , while the difference between them is the choice 
of the direction which separates the data into two subsets. 
By rotating the poles at I e [0° , 180°] and b e [-90°, 90°] 
in every 1° step, the original method uses definite num- 
ber of directions. With the developed method, one can 
have more random directions compared with the original 
one, by setting the repeating times (see below). There- 
fore we will take the developed method. The hemisphere 
comparison method involves the following steps. 

1. Generate a random direction 



(cos (f) sin 0, sin sin 9, cos ( 



(2.7) 



where </> G [0, 2tt) and 9 g [0, 27r] are random numbers 
with uniform probability distribution. Note that the uni- 
form probability distribution of these two variables (0, 9) 
will guarantee that every direction will be generated with 
the same probability, other priors on the distribution of 
(j) OT 9 may lead to a priori preferred direction. 

2. Split the dataset under consideration into two sub- 
sets according to the sign of the value fi-fidat, where hdat 
is the unit vector describing the direction of each SNIa 
in the dataset. Then we divide the data in two opposite 
hemispheres, denoted by up and down, respectively. 

3. Find the best-fitting values of {flmo,w) {(XlmO,wo) 
for CPL) on each hemisphere. Naturally, one can take 
use of the deceleration parameter to quantify the 
anisotropy level through the normalized difference 



Ago 
9o 



QO,u + lOA 



(2.8) 



where the subscripts u and d denote the up and down 
hemispheres, respectively. Note that the deceleration pa- 
rameter is a good diagnostic to quantify the anisotropy 
level. The hemisphere with larger go is expanding slower 
than the opposite. Larger normalized difference in one 
direction means that the anisotropy level in this direction 
is more notable. 

4. Repeat 400 times from step 1 to step 3, and find 
the maximum normalized difference defined by (2.8 1 for 



the Union2 data, thus one can obtain the corresponding 
direction of maximum anisotropy. 

As illustrated in [22], in order to maximize the effi- 
ciency, the number of directions should be no less than 
the number of data points of SNIa on each hemisphere. 
The reason for this is that changing the direction of an 
axis, does not change the corresponding 422. until a data 
point is crossed by the corresponding equator line. Such 
a crossing is expected to occur when the direction of an 
axis changes by the approximately mean angular separa- 
tion between data points. Thus, using more axes than 
the number of data points in a hemisphere does not im- 
prove the accuracy of the determination of the maximum 
anisotropy direction. Given that the number of data 
points per hemisphere for the Union2 dataset is about 



280, we have used 400 axes in our analysis, well above 
the value of 280. 



3. RESULTS 

Following the steps introduced in Sec. II, one can obtain 
the best-fitting values of go in each direction, and then 
find the maximum anisotropy direction. But in order to 
obtain the la errors of the maximal anisotropic direc- 
tion and the anisotropy level, we need to get the errors 
of the parameters ((fJmOjW;) foi' wCDM and {^IraOiW^) 
for CPL), which would propagate to qo and then to ^^a.^ 



because go is determined by w andTimO [see (2.3|] for 
wCDM model or Wq and fimo [see (2.5)] for CPL model. 



The analysis here is performed by using the Monte Carlo 
Markov Chain in the multidimensional parameter space 
to derive la errors on each hemispheres in the maxi- 
mum anisotropic direction. Accordingly, the la deviation 
from the maximum anisotropy level can be expressed as 
Ago _ Aqo^^aa: _|_ ^ corrcspondiugly the direction 
axes with la error can also be obtained. 

Further we explore the possible redshift dependence of 
the anisotropy. We implement a redshift tomography of 
the data and take the same procedure as before for all 
the following redshift sfices: 0-0.2, 0-0.4, 0-0.6, 0-0.8, 0- 
1.0. Our results for wCDM model are summarized in 
Table [T] And we also show the result for the full sample 
of Union2 (0 < z < 1.4) in Figure [l] 



redshift range 


I [degree] 


6 [degree] 


Ago 











- 0.2 1 332+f 1 -29;i2 




0-0.4 1 3i9;ii 1 -28;S 


1 oo^"-" 


0-0.6 1 30o;^3^ 1 -16;^:; 


Q4-"-ay 

'-'■^^+0.37 


0- 0.8 1 309;^' 1 2i;^^ 


91-"-^^ 

'-'■^^+0.30 


0-1.0 1 3ii;^^ 1 15;^? 


0-85+0.28 


0-1.4 1 314;^o' 1 28;?? 


^- ' ^+0.27 



Table 1: Directions of maximum anisotropy for several red- 
shift ranges of the Union2 data fitting with the wCDM model. 
The last column corresponds to the la errors propagated from 
Q.mO and w, which are obtained by MCMC method. 



The redshift tomography analysis here shows that 
the preferred axes are all located in a relatively small 
part of the North Galactic Hemisphere [around (Z, h) = 
(314°^2oo 1 28°^j^^o )], which is consistent with the result 
in |22j at la confidence level, thus indicates that un- 
der the assumption of wCDM model, the universe has 
a maximum acceleration direction. For different redshift 
slices, there are slight differences in the direction of pre- 
ferred axes, and the difference between the two oppo- 
site hemispheres is extremely obvious for the low-redshift 
slice {z < 0.2). 

Here we also give the best-fitting parameters of the 
wCDM model in Table [2j where the subscripts u and d 
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denote the up and down hemispheres. For the case with 
full SNIa data, see the last row in Table 2. Note that in 
the KCDM model, n^^^ = 0.30 and = 0.19 [22]. 



redshift range 








w'' 


0-0.2 


0.45 


0.10 


-0.80 


-0.51 


0-0.4 


0.44 


0.11 


-1.93 


-0.65 


0-0.6 


0.41 


0.11 


-1.91 


-0.70 


0-0.8 


0.37 


0.12 


-1.68 


-0.67 


0- 1.0 


0.37 


0.18 


-1.68 


-0.76 


0- 1.4 


0.35 


0.22 


-1.46 


-0.77 



Table 2: Best-fitting parameters of the wCDM model for 
different redshift slices on the opposite hemispheres in the 
direction of the maximum anisotropy, where u denotes the 
hemisphere corresponding to larger accelerations, while d de- 
notes the opposite hemisphere. 

Note that 17^0 is used as the diagnostic to the 
anisotropy level in [22] . Here we use qo as the diagnostic 
to the anisotropy level. We found here that qo is more 
sensitive to the errors of the data, since go is relative to 
the second order derivative of the luminosity distance. 
And also because more parameters involved in the anal- 
ysis will introduce more uncertainties when we try to de- 
termine the error of qo , therefore the errors here are larger 
than those in [22]. For example, the dots colored accord- 
ing to the sign and magnitude of the anisotropy level will 
be more scattered on the unit sphere, as shown in Fig- 
ure [1} But the meaning of using qo is obvious that the 
hemisphere with larger go is expanding slower than the 
opposite, and that the direction with normalized differ- 
ence ^|2o jneans that the anisotropy level in this direction 
is more notable. What's more, it is clear from Figure [T] 
that the sphere is divided into two distinct hemispheres, 
one with smaller accelerations and the other with larger 
accelerations. 

With the same procedure we also study this issue 
for a dynamical dark energy model with the CPL 
parametrization. By use of all Union2 data points, 
we find that the direction of preferred axis is (/, b) = 
(309°^30o , 21°;jl||o ), and the corresponding maximum 
anisotropy is = 0.76^o;4^. The result shows that 

it is not much different from the case of the wCDM 
model. This means that the best-fitting values of the 
preferred direction is not very sensitive to the dark en- 
ergy models. 



4. CONCLUSIONS 

From some astronomical observations and some the- 
oretical models of the universe, there seemingly exists 
some evidence for a cosmological preferred axis jlOj . If 



such a cosmological preferred axis indeed exists, one has 

Aq/q„=0.79 




Aq/q„= -0.79 

Fig. 1: The dots on the unit sphere colored according to the 
sign and magnitude of the anisotropy level are the directions 
of the random axes. The hemisphere shown on the left panel 
is the one corresponding to smaller accelerations with the 
preferred axis {l,b) = (134°, -28°) [qo = -0.40 in this di- 
rection], while the right corresponds to the one with larger 
accelerations and the preferred axis {l,b) — (314°, 28°) and 
qo = -0.92. 



to seriously consider an anisotropic cosmological model 
as a realistic model, instead of the FRW universe model. 

In this paper we investigated the existence of 
anisotropy of the universe by employing the hemisphere 
comparison method and the Union2 SNIa dataset and 
found this preferred direction. We used the present 
value of the deceleration parameter go to quantify 
the anisotropy level of the two hemispheres. For 
the wCDM model, the preferred direction is (Z,6) — 
{?>lA°-^.,2i°-li:), and ^ = O.Tg;^:^^. While in 
the case of CPL model, the direction of preferred axis 
is (Z,5) = (309°^|o° ' 21°^|5o ), and correspondingly the 
maximum anisotropy level is = 0.76^o4^. Com- 

paring with the result given in [22], where the KCDM 
model is employed, our results are basically in agreement. 
This means that the best-fitting preferred direction is not 
much sensitive to the dark energy models. 

Finally let us notice that it can be seen from Table 
1 that the result is weakly dependent on redshift if the 
redshift tomography analysis is employed. 
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